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We theoretically investigate effects of quantum fluctua- 
tion between an even and an odd charge state on the trans- 
port properties of a normal-superconducting-normal single- 
electron tunneling transistor. The charge fluctuation is dis- 
cussed beyond the orthodox theory. We find that the charg- 
ing energy renormalization enhances the parity effect: The 
Coulomb blockade regime for the odd state is reduced and 
that for the even state is widened. We show that the renor- 
malization factor can be obtained experimentally and that the 
renormalization effect is weakened by applied bias voltage. 



I. INTRODUCTION 

Quantum fluctuations play important roles on trans- 
port properties of mesoscopic systems at low tempera- 
tures. Recently the charge fluctuation in Coulomb is- 
lands has attracted much attention. There have been 
much development in the strong tunneling regime. A 
theoretical prediction, a rcnornwlization of the conduc- 
tance and charging energy £yjT3, has been confirmed 
experimentallyBO. 

When the island is made of superconductor, the 
"parity effect" plays an important role on transport 
propertiesETEI. When Ec exceeds the superconducting 
gap A, the parity effect appears in the period of Coulomb 
oscillation where intervals are either elongated or short- 
ened for even or odd occupancy of electrons respectively. 
At a resonance, an even and an odd state are degenerate 
and an unpaired "odd" electron makes dominant contri- 
bution to the tunneling current. For Ec < A, odd states 
are no longer stable for every gate voltage and thus 2e- 
periodic Coulomb oscillation appearsEl 

Therefore, it is an intriguing question to ask how the 
parity effect with charge fluctuation affect on transport 
properties. In this paper, we investigate the conduction 
properties of a superconducting island for Ec > A with 
a special attention to the charge fluctuation between an 
even and an odd state. For simplicity we limit ourselves 
to the discussion at zero temperature limit aad neglect 
the interference effect around tunnel barriersta. We find 
that the interval of conductance peaks related to the 
Coulomb blockade (CB) regime for the odd state is short- 
ened and that related to the even state is widened as a 
result of the charging energy renormalization. We show 
that the renormalization factor can be obtained experi- 
mentally and that the renormalization effect is weakened 
by applied bias voltage. 

The outline of this paper is as follows. In Sec. ||, 
we introduce the model and derive the generating func- 



tional in the path-integral representation. We als o pro - 
pose the approximate generating functional Eq. (2.12). 
In Sec. Ill, we derive the current expression by using the 
functional derivative. In Sec. IV we show numerically 
evaluated results and have some discussions. Section |v| 
summarizes our results. 



II. MODEL AND GENERATING FUNCTIONAL 

Figure |l](a) shows an equivalent circuit of a normal- 
superconducting-normal (NSN) transistor. A supercon- 
ducting island exchanges quasiparticles (QPs) with a left 
(right) lead via a small tunnel junction characterized by 
the tunnel matrix element Tl(r) an d a capacitor Ctyr) 
and is coupled to a gate via a capacitor Cq. In the follow- 
ing discussion, we limit ourselves to the symmetric case, 
Cl = Cr and Tl = Tr. We use the two-state model 
to describe the strong Coulomb interaction. We assume 
that there are even (odd) number of electrons in a charge 
state |0(1)>. 
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FIG. 1. (a) The equivalent circuit of the NSN transistor, 
(b) The normalized density of states in the superconducting 
island. 

The total Hamiltonian H is given by the sum of Hq 
and a tunneling Hamiltonian Ht, which is adiabatically 
switched on at remote past and off at distant future. Hq 
is defined as 



Hq — ^ £r/c a lkn a rkn + A 



<J Z + 1 



E c (Q G /ef 



r.k.n 



(2.1) 



where a r fe„ is the annihilation operator of a QP in the 
lead (r = L, R) or in the island (r = I) with wave vector 
k and transverse channel n which includes spins, a is 
t he effective spin-1/2 operator. The second term of Eq. 
fl2.l| ) describes the energy difference between two charge 
states. The excitation energy Ao = Ec(l — 2Qq/c) + A 
i s re lated to a gate charge Qg- The third term of Eq. 
(|2.l|) represents the charging energy for the even state. 
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The density of states (DOS) for normal metal leads is 
approximated as constant: p r (e) = ^2 k S(e — £ i fe) = P* 
(r = L, R). Figure |l|(b) shows the normalized DOS in 
the island pi(e)/ where pj is DOS for the normal state. 
We assume that an unpaired electron occupies the lowest 
energy level above the superconducting gap in the limit 
of K: 
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(2.2) 



where we consider the condition that the charge and the 
spin relaxation in the island are enough fast and occupa- 
tion probabilities for the up- and the down-spin unpaired 
electron are both 1/2. We expect that the most dominant 
part of the superconducting correlation effects can be in- 
cluded via the shape of DOS, because any other process 
than the QP tunneling, such as the Andreev reflection 
process, requires higher energy ~ Ec- 
The tunneling Hamiltonian 



An 



E E * 

r=L,Rfc,fc',n 
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a Ikn a rk'n<J+ 



h.c. 



(2.3) 



describes the electron tunneling across the junctions and 
simultaneous flip of the effective spin. The phase dif- 
ference between the lead r and the island is written as 
tp T = K v ip. ip is the phase difference between two leads 
and is regarded as a source field. kl(r) = 1/2 (—1/2) 
characterizes the voltage drop between the left (right) 
lead and the island. 

The Hamiltonian includes the effective spin-1/2 opera- 
tor a and thus the Wick's theorem for fermions or bosons 
cannot be used. To circumvent this.dijawback, we employ 
the drone-fermion representationli3o, which is the map 
of g onto two fermion operators c and d as <r+ = c^(j> and 
<j z = 2c)c — 1, where <fi = eft + d is a Majoranapfermion 
operator. In our calculation, we reformulate Ref.EJ based 
on the Schwinger-Keldysh approacniaij which enables us 
to obtain the average current, the current noise and any 
higher order moments systematically by the functional 
derivative m tprms of the phase difference, i.e. a vector 
potentialBEa. 

The action can be obtained from the total Hamiltonian 
H by following the standard procedure^: 



S [c ,c,d , d, a lkn , a lkn ] 

= / dt{c(t)*(ihd t - A )c(t) + ihd(t)*d t d(t) 
Jc 

+ ^ a rkn{t)*(iftd t - e lk )a lkn (t) 

r,k,n 

+ T r e l ^a rkn {tya lkln {t)v+(t)+\i.c.}, 



(2.4) 



r=L,R 
k,k'n 



where the derivative and the integration are performed 
along the closed time-path C consisting of the forward 



branch C+, the.backward branch C_, and the imaginary 
time path C* r El(Fig. |). The dc grees of freedom for ip 
are duplicated, i.e., we can define (p + and ip~ on the 
forward and the backward branch, respectively. a lkn , c 
and d are Grassmann variables which satisfy the anti- 
periodic boundary condition. By tracing out QP degrees 
of freedonJa (Appendix. |A|) , we obtain the following ef- 
fective action for the c- and d-field: 

S [c* , c, d* , d] 



= / dt{c(t)*(ihd t - A )c(t) + ihd(t)*d t d(t)} 
Jc 

+ [ dt 1 dh{a + (t 1 )a(t 1 ,t 2 )a-(t 2 ) + 0(T r 4 )}, 
Jc 



(2.5) 



where trivial constants are omitted. The second inte- 



gral of Eq. (2.5) describes the tunneling process, a — 
S r =L r a r is a particle- hole Green function (GF), written 
as 



a T {t,t') 
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(2.6) 
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Here N c h is the number of transverse channels and g lk 
is a free QP GF in the lead r (r = L, R) or in the island 
(r = I). The inverse of g rk is defined as 



97k M = (ihd t -s lk )S(t,t'), 



(2.7) 



where ^-function is defined on C and g Tk satisfies the 
anti-periodic boundary condition: g lk (t,—oa € C+) = 
-g r k(t,-ihp - oo). 

In the wide junction limit, A^ c h — > oo, the terms higher 
than Tj 4 , which describes elastic co-tunneling processes, 
can be neglectedu. The particle-hole GF a describes tun- 
neling and relaxation process related to the lowest order 
sequential tunneling. However, by tracing out c- or d- 
field, a number of a are coupled. Therefore, the effective 
action describes the higher order inelastic co-tunneling 
processes, too. 

Next we trace out drone-fermion fields. Firstly, we 
trace out the degrees of freedom for c- field. The result- 
ing effective action for d-field includes the many-body 
interaction, and cannot be solved exactly. Thus we in- 
troduce a liner source term J c dtiJ(ti)4>(ti), where J is a 
Grassmann variable defined on C. By tracing out d-field 
degrees of freedom, we obtain the generating functional 
as 



Z = exp - J2 



-Tr 



9c Sl a 6J 



x exp 



2ih 



dt 1 dt 2 J{ti)g 4> (ti,t 2 )J(t 2 , 



,7=0 

(2.8) 



where trace is performed over C and products in the trace 
represent the integration along C. We omit the factor 2 



2 



which is the partition function of d-field. Here the c-field 
and the d-field GFs are defined as 



9c' (*,*') 



(ihd t -Ao)S(t,t'), 
ihd t S{t,t')/2. 



(2.9) 
(2.10) 



Both GFs satisfy the anti-periodic boundary condition. 

The generating functional for connected GF, W — 
—ihlnZ, is evaluated by performing the perturbation 
series expansion in powers of a, namely the dimen- 
sionless junction conductance ao — X) r =L r a r wnere 
a° = iV c hT r 2 pjp r . For example, the first order contri- 
bution is written as W^ 1 = ihTr[g c Y< c ], where the self- 
energy S c = 53 r=L R E r is defined as 



S r (t,t / ) = —ihg^t' ,t) a T (t,t'). 



(2.11) 



A finite order contribution causes a divergence of the 
physical quantity, such as the average charge, at the de- 
generacy point Ao = (Appendix. |B|). To regularize the 
divergence, one must take infinite orders into account. 
The most simple way is to sum up order-n c-field correc- 
tions (<7 C EJ^\ This strategy is also adopted by Isawa et. 
al. in Ref.E3, though their formulation is different from 
ours. The resulting approximate generating functional is 
expressed as 



V7 = -ifiTr[lnG- 1 ], 



(2.12) 



where G c is the full c-field GF whose inverse is defined 
as 



G- 1 (t,t')=g- 1 (t,t')-i: c (t,t'). 



(2.13) 



As we show in the next section, the approximate gen- 
crating functional formally reproduces the-result of the 
resonant tunneling approximation (RTA)0II3. 
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FIG. 2. The closed time-path going from — oo to oo (C+), 
going back to — oo (CL), connecting the imaginary time path 
CV and closing at t — —oo — ih/3. 



III. AVERAGE CURRENT 

The tunneling current is obtained by functional deriva- 
tive of the generating functional W with respect to 
the phase differencetZl as I = {e/K)5W / 8ip A (f) \ v& _ 

Ei=lr k '-'i' The relative coordinate — ip + — tp_ is 
a fictitious variable and must be 0. The center-of-mass 
coordinate, (p c — (y>_|_+y_)/2, is a physical variable fixed 
at eVt/h. By regarding (p T as the formally independent 
variable, the tunneling current through the junction r ,/ r , 
is expressed as 



7 r (t) = 



SW 



-e Tr 



5(p T 



G, 



VA=0 



G c E r 



5 (fir 



(3.1) 



Here, we pay attention to the fact that only the self- 
energy depends on the phase difference. Next we project 
the time defined along C to the real axis. As the tunnel- 
ing Hamiltonian is turned on and off adiabatically, the 
particle-hol e G F is zero on the imaginary time path C r . 
Thus, Eq. Q is rewritten as, 



eTr 



G c T 1 I] r 



S(fi T 



G r 



Stp lA (t) Sip rA (t) 



-E.T 1 



, (3.2) 



where the trace is performed in the 2x2 Keldysh space 
and products represent the integration along the real 
axis. x s (s = 0, 1, 2, 3) is the Pauli matrix in the Keldysh 
space. Here GF and the phase difference denoted with 
tilde are the 2x2 matrix for GF and -that for the scalar 
variable in the physical representatior&B: 



G <: 



<Pr = Q 



G^ 
r<R r<K 



ip T+ 








f 



<PtA- 



(3.3) 



(3.4) 



where GFs with superscripts, A, R and K, represent the 
retarded, the advanced and the Keldysh component, re- 
spectively. The practical calculations of these compo- 
nents are shown in -Appendix. The matrix Q is 
the Keldysh rotatort£l. From Eq. ( |3.4[) , we derive the 
useful relation for the functional derivative technique: 
5(p. v {t')/8ip vA (t) = T°5{t - t')/2. By using the prop- 
erty of GF in the physical representation G(t,t'y = 



5 G(i',i)r 3 



T T 



we can see that the second term of Eq. (3.2) is minus the 
complex conjugate of the firs t te rm. By performing the 
Fourier transformation, Eq. (3.2) becomes as 



2eRe / cftiTr 



-01 



= 2< Re / ylv 



GefttijT^fa,*) — 



(3.5) 



where G? = G? 



Gg , etc 



By usin g the expressions for 
the full c-field GF (Eqs QCl^ ) an d (jci7| )) and those for 



the self-energy (Eqs. (C14) and (C15)), we obtain the 
expression for the average current in the limit of K: 



I = 



Gk 

e 



:V/2 



de 



eV/2 



(c) 



2iImGf(e), (3.6) 
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where Gk = e 2 /h is the conductance quantum and af 
is the Keldysh component of the particle-hole GF (Eq. 
(|C9|)). The particle-hole GF of the superconducting is- 
land in Eq. (il) is given by af(e) — — 27ri a° \p(5e r — A)| 
where Se r — e — n r eV. The spectral function p(e) is given 
by 



-Ve 2 
p(e) = { l/2p° 



A 2 



l/2p° + 



(-£> < e < - 
(\e\ < A) 
(A < e < D) 



A) 



(3.7) 



where ther-,high energy cut-off D (= Ec S> A) is 
introducedcl. The imaginary part of the retarded c-field 
GF Gf (e) = l/(e- A - Sf (e)) describes the excitation 
spectral density of the charge state. The self-energy of 
c-field is expressed by (e) = J2i=l r a rR{5£ T — A) — 
17(e), where the function R is written as 



R(e) 



x Re 



2 e In (2D/A) + 2Ve 2 - A 2 
Be 



tanh 



A 2 - VD 2 - AVe 2 - A 2 



ln(|( e -D)/(e + A)\)/p° 1 , 



He) 
(3, 



where h(e) = for |e| < A and h(e) = tantT 1 (£>/e) 
otherwise. The imaginary part of the self-energy 7(e) = 
— Im \a (e)] /2 represents the life-time broadening of the 
charge state caused by dissipative charge fluctuation. 
The charge fluctuation described by a K is suppressed 
in the low energy range, < e < 2A <C D, because for 
the odd state, there are no states for QPs in the range 
— 2A < e < as shown in Fig. 0(b). The asymmetry 
in particle-hole GF a K (e) causes the renormalization of 
the peak position of ImGf (e). ft should be noted that 
the origin of the renormalization is different from that 
for the normal metal islandBTJ and the double-islandEll 
where particle-hol e GF is symmetric. 

The expression ( |3^6| ) is formally equivalent to that ob- 
tained within RTAca The validity oLRTA, which is de- 
veloped for the normal metal islandou, is not obvious 
when the superconducting correlation exists. However, 
we expect that RTA can be used as a start point of ap- 
proximation when QP tunneling is the main transport 
mechanism. The expectation is partly justified by the 
following two results. First, pEq. ( p.6| ) reproduces the 
result of the orthodox theoryEil in the limit of small di- 
mensionless junction conduct ance ao- Secondly, in CB 
regime (\eV\ <§; |A |), Eq. ( |3.6| ) is approximately ex- 
pressed as 

G K (2^) 2 <4(2eA^) 1 / 2 T//3 P °A 2 , 

which is equal to the. previous expression obtained by 
Averin and NazarovEa in CB regime for the odd state 
(\eV\ — Ao) except for the renormalization of the peak 
position 



An = An 



2a Aln(2£; c /A) - ln^p^c /ttq ) / pi 



As for the CB regime for the even state i\eV\ <C Ao), our 
result does not reproduces the previous result that the 
inelastic co-tunneling current is zeroEil. Here we should 
comment on the limit of the application of RTA. RTA 
takes account of high-order inelastic co-tunneling pro- 
cesses where at most one particl&.hole excitation is cre- 
ated inside the island at a moment! 2 ]. And thus, the simple 
application of RTA for the even state, where a low-energy 
particle-hole excitation cannot be created inside the is- 
land, causes us to count some undesirable processes. To 
improve the approximation in CB regime for the even 
state is outside the scope of this paper. 



IV. RESULTS AND DISCUSSIONS 

In our numerical calculation, we choose the supercon- 
ducting gap and DOS as A/E c = 1CT 1 and p^E c = 10 3 , 
respectively. The variable p® A = 10 2 is the same order 
as that of the Al island whose volume is 10 4 nm 3 . We con- 
firm that our numerical calculation reproduces the spec- 
tral sum rule of the full c-field GF J delmG^(e) = — it 
within 0.2% accuracy. 
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An/Ec 

FIG. 3. The excitation energy dependence of differential 
conductance for an = 0.1 (solid line), 0.2 (dashed line) 0.3, 
(dot-dashed line) and 0.4 (dotted line) at eV/E c = 10~ 4 . 
Inset: A schematic diagram of the total energy with (dashed 
lines) and without charge fluctuation (solid lines). 

Let us first consider the effect of the charge fluctuation 
on the conductance G. Using Eq. ( |3.6| ), we calculate G 

as 

G/Go ~ a (l/2p°) 2 /{A 2 + (na Q /2fi ) 2 }, 

where Go = Gk(^) 2 <XqOq /ao- Figure || shows the exci- 
tation energy dependence of differential conductance at 
small bias voltage (eV/Ec = 10~ 4 ) for various ao- The 
width of the conductance peak is approximately given 
by 7rao I^Pi inversely proportional to the life-time of the 
unpaired electron. The conductance peak shifts leftward 
as ao increases. The peak position, namely the degen- 
eracy point, is at A ~ — 2aoAln(2i?c , /A). At the de- 
generacy point, the gate charge Qc/e is approximately 
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A/{2zE c ) + 1/2, where z = 1/(1 + 2a ln(2£' c /A)) is 
the renormalization factor. The charge fluctuation re- 
duces the charging energy, while it does not affect on 
the magnitude of the superconducting gap. This result 
is qualitatively different from that of the normal island, 
where both of the charging energy and the conductance 
are renormalizedel. 

The origin of the shift of the degeneracy point can be 
understood by looking at the schematic diagram of gate 
charge dependence of the total energy shown in the inset 
of Fig. [| The dashed and the solid line indicate the total 
energy with and without charge fluctuation, respectively. 
Here we assume that the renormalization factor is inde- 
pendent of Ao. One can see that the curvature of the 
energy curve with the charge fluctuation is smaller than 
that without the charge fluctuation since the charging 
energy is reduced due to the charge fluctuation. Thus 
the degeneracy points are shifted to reduce CB regime 
for the odd state. 

It should be stressed that the renormalization factor 
can be determined experimentally. From the Coulomb 
oscillation, one can obtain the ratio of the superconduct- 
ing gap to the renormalized charging energy A/ (zEq) = 

(Qq — Qq) I {Qq + Qq) , where Qq^ is the interval corre- 
sponding to odd (even) state. Here the "bare" charging 
energy Ec and the superconducting gap A can also be 
obtained experimentally. For example, one can obtain 
Ec from the temperature dependence of the conductance 
peak for the normal stateS. The superconducting gap A 
can be obtained from the I-V characteristic as discussed 
below. Therefore, one can estimate the renormalization 
factor quantitatively by analyzing the experimental re- 
sults. 
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FIG. 4. The excitation energy and the bias voltage depen- 
dence of the tunneling current for (a) ao = 0.1 and (b) 0.3. 
The corresponding density plots of differential conductance 
are shown in panels (c) and (d). 



Figures ||(a) and [|(b) show the excitation and bias 
voltage dependence of the tunneling current for ao = 0.1 
and 0.3, respectively. The corresponding density plots 
of differential conductance are shown in Figs. [|(c) and 
^(d). Dotted lines indicate boundaries of four regions 
predicted by the orthodox theoryEx (1) CB regime for 
the odd state, (2) the unpaired electron tunneling regime, 
(3) CB regime for the even state, and (4) the QP tunnel- 
ing regime. The tunneling current shows the plateau in 
the region (2). As a becomes large (Figs. |](a) and||(b)), 
one can see that the structure of the tunneling current 
is smeared by the life-time broadening, especially around 
the QP tunneling regime (4) . However, the threshold bias 
voltage where a QP tunneling channel opens eV = 2A 
does not depend on ao and thus one can obtain A exper- 
imentally. 

The striking feature in Figs. |(c) and |(d) is that at the 
boundary of CB regime for the even state, the conduc- 
tance peak survives until eV reaches 2A. This is because 
the dissipative charge fluctuation is blocked due to the 
superconducting gap and only the unpaired electron con- 
tributes to the life-time broadening: 7 ~ ira /2p®. On 
the other hand, at the boundary of CB regime for the 
odd state, the conductance peak is gradually smeared. 
This is because the dissipative charge fluctuation in- 
creases with increasing bias voltage, and the life-time 
broadening 7 is proportional to the number of QP states: 
7 ~ iraoV A eV. 

The renormalization effect is expected to be weakened 
by applied bias voltage because eV gives the low energy 
cut-off for the renormalization factora. The reduction of 
the renormalization effect can be deduced by measuring 
the peak position of the conductance at the boundary of 
CB regime for the even state (Figs, ^(c) or ||(d)). In 
the limit of zero bias voltage, the position of the conduc- 
tance peak deviates from that predict by the orthodox 
theory(dotted line). As the bias voltage increases, the 
position of the peak approaches the dotted line and they 
meet at eV = 2A. It means that the renormalization 
effect is weakened with increasing bias voltage. 

V. SUMMARY 

In conclusion, we have theoretically investigated the 
transport properties of NSN transistor showing the parity 
effect and the quantum fluctuation of charge. We con- 
sidered the quantum fluctuation between the even and 
the odd state caused by the inelastic resonant tunnel- 
ing precess. We found that the charge fluctuation causes 
the renormalization of the charging energy, and that CB 
regime for the odd state is reduced. The renormalization 
factor can be obtained experimentally. The renormal- 
ization effect is weakened by applied bias voltage. At 
the boundary of CB regime for the even state, the con- 
ductance peak is robust against the charge fluctuation 
because only one unpaired electron can contribute to the 
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life-time broadening of the charge state. On the con- 
trary, at the boundary of CB regime for the odd state, 
the conductance peak is smeared since the dissipative 
charge fluctuation increases with increasing applied bias 
voltage. 



where / + (e) is defined with the Fermi function f~(e) = 
l/(e /3e + 1) as / + (e) = /~(e)e^ e . In the opposite case, 
t < t', g(t,t') is written as 

-f-(s)e X p(e(t-t , )/(ih))/(ih). (A4) 
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APPENDIX A: FUNCTIONAL INTEGRAL ON 
THE CLOSED TIME-PATH 

In this Appendix, we demonstrate the method to cal- 
culate the functional integral defined on the closed time 
path C. For generality, we consider the following action 



S 



dt{c{t)*{ihd t - e)c{t) + J{t)*c{t) + h.c.}, (Al) 



where c is a Grassmann variable. A complex variable 
J is a source field and is zero on C,. In order to uti- 
lize the time-slicing technique on Cc3, we introduce the 
cut-off time to > which restricts the range of C± 
from — to to to and we take the limit as to — > 00 a f~ 
ter all calculations. We divide contours C+, C_ and 
C T into N pieces, respectively. The discretized time is 
f « = Z)i=i, n e * ~ *o where e„ = 2t /N for n = 1, • • • , N, 
e n = -2t '/N for n = N + 1, • • • ,2N and e„ = -ih/3/N 
for n = 2N + 1, • • • , 3 AT. The action is written in the 
discretized form as J2^!j=i c *^y' c j + !Ci=i e iJ* c i + h.c, 
where Ay = ihSij — (ih + eie)6' i j +1 . 5[ - +1 is the Kro- 



for 



necker's delta Si j+i for j = 1, • • • , 3N — 1 and — Si i 
j = 3N. The integration is performed in the same way as 
the imaginary time path-integralES. By taking the limit 
as to — > oo after taking the continuous limit N — > oo, we 
obtain the generating functional as 



Z = Zo exp 



dtdt'J*(t)g(t,t')J(t') 



c 



(A2) 



Here Zo — expTrln(g x ) becomes the partition func- 
tion for fermion, 1 + e~ ,3e , because the contribution from 
C+ and that from C_ chancel each other and only the 
contribution from C T r ema ins. 

Though GF in Eq. ( |A2p is defined on C, in the prac- 
tical calculations GF defined on C+ + C_ is needed. If 
t > t' with respect to the contour C + + C_, g(t,t') is 
written with the time projected on the real axis as 



f+(e)exp(e(t-t')/(ih))/(iH), 



(A3) 



APPENDIX B: FIRST ORDER PERTURBATION 
THEORY 

In this Appendix, we demonstrate that the lowest order 
perturbation theory gives the divergent average charge 
at the degeneracy point within our formulation. For 
simplicity, we consider the case of the equilibrium state 
and in the limit of zero temperature. In this case, the 
imaginary-time formalismE£l is convenient for calculations 
and the grand canonical potential is evaluated perturba- 
tively in the same way as Sec. ||. The first order term 
of the grand canonical potential f^ 1 ) becomes the similar 
form as and consists of the thermal GF for c-field, 
field and the particle- hole given by l/(«w; — A ), 2/(iwi] 
and a(iv n ), respectively. Here, u>i and v n are the fermion 
and the boson Matsubara frequency, respectively. 



i 



(3 2 f-' iv n + iui - A iuJi 



a(iu n ) 



-ilm[a(e-M0)l 

An — £ 



where P denotes the Cauchy's principal value integral 
and N~(s) — l/(e^ £ — 1) is the Bose distribution func- 
tion. Here the numerator of the integrand is the spectral 
function of the particle- hole GF written as aop{e — A) 
where p is given by Eq. (3/7). The correction of the aver- 
age charge is evaluated by the derivative of Cl^ in terms 
of the excitation energyo. In the limit of zero tempera- 
ture, it becomes as dQ^/dAo ~ iop_sgn(A ) dR(A — 
A)/dA where R is given by Eq. (|3.8|). As a result, in 
the CB regime for even state < An. <C Ec, we can 
see that the correction diverges as ~ ceo^y/ A/ (2Ao)/2, 
wher e we omit the contribution form the last term of Eq. 
( |3.8[ ), which is negligible when we consider the life-time 
effect. For normal state, our formulation reproduces the 
well known result, the log-divergenceLll The divergence 
for NSN transistor is strong as compared with the log- 
divergence. 



APPENDIX C: GREEN FUNCTIONS 

In this Appendix we calculate GFs. Before proceed- 
ing to the calculation of each GF, we demonstrate the 
method to obtain GF by solving the differential equation 

g- 1 (t,t') = (iHd t -e)S(t,t'), (CI) 
g{t, -to G C+) - -g(t, -ihf3 - to), (C2) 



G 



as compjj 
pendix 
pendix 
t, t' e Cr 
the relation 



entary to the method demonstrated in Ap- 
Here the definition of to is given in Ap- 
First, we calculate GF defined on C T , viz. 
In this case, the ^-function formally satisfies 

-tft/9 



J " >IJ dt5(t) = 1. By using the anti-periodic 
boundary condition Eq. (|C2|), we-can solve Eq. (CI) in 
the same way as the thermal GFc3. From the solution, 
we can show the following relation as 



?(-to G C T , -t Q G C±) = ±f ± (e)/(ih), 



(C3) 



where, we use conditions g(—to e C_,— to G C+) = 
-fli(-*o G C_, -t - ifi/3) and flr(— *q G C+, -t G C*_) = 
—g(—to — ifi/3,—to& fy)- Second, we calculate GFs de- 
fined on the real axisea. By projecting the time defined 
on C±, onto the real axis, g(t, t') is projected as g ±=t (t, t') 
for t,t' G C± and as g ±T (t,t') for f G C± and t' E C T . 
It should be careful that the arguments £ and t' of GFs 
with superscripts are the time projected onto the real 
axis. These four GFs satisfy the following differential 
equations as 



(ihd t - e) 



g^W) = ±6(t-t>) 
g ±T (t,t>) = o, 



where the ^-function is defined on the real axis. These 
differential equations can be solved by using Eq. (C3) 
and matching conditions, g s± (t,to) = g ST (t,to) and 
g ±s (to,t') = g TS (t ,t') where s = + or — . The resulting 
GFs are 



g++(t,t') = g-+(t,t')9(t-t')+g+-(t,t')0(t'-t) 
g—(t,t') = g-+(t,t')9(t' -t)+g+-(t,t')9(t-t') 



(t,f) = ±f ± (e)exp(e(t-t>)/(ih))/(ih). 



The results are equivalent to Eqs. (A3) and (A4). 

The four GFs are-aomponents of 2 x 2 GF in the single 
time representation^. It is known that this representa- 
tion includes the redundancy which can be removed by 
the Keldysh rotatiorJlaij. After the Keldysh rotation, 
we obtain 2x2 GF in the physical representation (see 
Eq. ( |3.3| )). In the following discussions, we summerize 
the retarded and the Keldysh components of 2 x 2 GF in 
the physical representation. The advanced component is 
obtained by taking the complex conjugate of the retarded 
component in the energy space. 

The differential equation Eq. (|]7]) can be solved in the 
same way as above example: 



l/(e + irj- e Tk 
\2T 



-2wrtanh | 



) 6(e-e lk ), 



(C4) 
(C5) 



where, 77 is a positive small value and the 5-function in 
the energy space is defined as 6(e) 
By using Eqs. and 



7T € 2 +7} 2 

we can calculate the 



retarded and the Keldysh components of particle-hole GF 
defined by Eq. (|j]) . In the following, we fix the relative 
coordinate of the phase difference as <pa(£) = 0, and the 



center-of-mass coordinate as <p c (t) — eVt/H.j- The loop 
diagram can be calculate in the standard wajc3: 



af(s)=iV ch T r 2 ^ f | 

u 1.1 J 



k,k' 



de[_ 

2l7T 



= -iira°p(6e r - A), 
*?(e)=N ch T?Y: IP- 

h hi J 



k,k' 



gf fc (fe r + e') gi K k ,(e') - <4(fe r + e')g^(e') 
2 

-2*™?p(fe r -A)coth(^ , 



(C6) 
(C7) 



(C8) 
(C9) 



where the spectral function p(e) is defin ed b y Eq. (3.7). 

By solving the differential equations ( |2.9| ) and ( 2.10 ) 
we obtain GFs for c- and those for <i-field: 



g$(e)=2/(s + ir ] ), 
4(e) = 0, 

9 f(e) = l/(e + i 7 ? -A ), 



f (e) = -2*7rtanh(— j«5(e-A ). 



(CIO) 
(Cll) 
(C12) 
(C13) 



The retarded and the Keldysh component of self- 



energy for c- field Eq. ( 2.11 ) are calculated by using Eqs. 
©, ©, O, and IClll): 



(C14) 
(C15) 



^R( c \ - ( de ' ia ?( £/ ) 

^(e)=a?(s)-a?(e) = 2a?(sl 

whe re we utilize the similar expression as Eqs. ( p6| ) and 
©• 

The expressions for full c-field GF i s obt ained by solv- 
ing the inverse Dyson equation Eq. (2.13). By project- 
ing the time on C± onto the real axis, and performing 
the Fourier transformation, we obtain the matrix Dyson 
equation as G c (e) — g c (s) — 3c(£)T 1 S c (e)r 1 G c (£). The 
matrix Dyson equation is solved easily: 



G c fl ( £ ) = l/( £ + ^-A -Ef( £ )) ) 



(C16) 



Gf (e) = G* (e) {sf (e) - 2* r^tanh (JL) } G*(e), (C17) 

where we use the definition of <5-function in the energy 
space. Taking the limit 77 — > 0, we obtain the final form 
for full c-field GF. 
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